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We study chiral symmetry breaking in quenched strong-coupling QED4 in arbitrary covariant gauge within 
the Dyson-Schwinger equation formalism. A recently developed numerical renormalization program is fully im- 
plemented. Results are compared for three different fermion-photon proper vertex Ansatze: bare 7 M , minimal 
Ball-Chiu, and Curtis-Pennington. The procedure is straightforward to implement and numerically stable. We 
discuss the chiral limit and observe that in this limit the renormalized axial current is conserved. A detailed study 
of residual gauge dependence due to the vertex choice is in progress. The relevance for lattice studies is discussed. 
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1. INTRODUCTION 

Strong coupling QED in three space and one 
time dimension has been studied within the 
Dyson-Schwinger Equation (DSE) formalism for 
some time [fj]J^,|3| . For a recent review of Dyson- 
Schwinger equations and their application see for 
example Ref. [§. DCSB occurs when the fermion 
propagator develops a nonzero scalar self-energy 
in the absence of an explicit chiral symmetry 
breaking (ECSB) fermion mass. We refer to cou- 
pling constants strong enough to induce DCSB as 
supercritical and those weaker arc called subcrit- 
ical. We write the fermion propagator as 



w tf-M(p 2 ) A{p 2 ) i>- B{p 2 ) w 

with Z(p 2 ) the finite momentum-dependent 
fermion renormalization, and B(p 2 ) the scalar 
self-energy. In the absence of an ECSB bare 
electron mass, by definition DCSB occurs when 
B(p 2 ) ^ 0. Note that A{p 2 ) = 1/Z{p 2 ) and 
M{p 2 ) = B(p 2 )/A(p 2 ). 

Many studies, even until quite recently, have 
used the bare vertex as an Ansatz for the one- 
particle irreducible (1-PI) vertex Y v (k : p) despite 
the fact that this violates the Ward-Takahashi 
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Identity (WTI) fj. With an Y 01 tnese Ansatze 
the resulting fermion propagator is not gauge- 
covariant, i.e., physical quantities such as the crit- 
ical coupling for dynamical symmetry breaking, 
or the mass itself, are gauge-dependent |J. A 
general form for T"(k,p) which does satisfy the 
Ward Identity was given by Ball and Chiu in 
1980 0; it consists of a minimal longitudinally 
constrained term which satisfies the WTI, and a 
set of tensors spanning the subspace transverse 
to the photon momentum q. Although the WTI 
is necessary for gauge-invariance, it is not a suffi- 
cient condition; further, with many of these ver- 
tex Ansatze the fermion propagator DSE is not 
multiplicatively renormalizable. There has been 
much recent research on the use of the trans- 
verse parts of the vertex to ensure both gauge- 
covariant and multiplicatively renormalizable so- 
lutions §,§|l^@||,yi, some of which will be 
discussed below. 

What was common to essentially all of the pre- 
vious studies is that the fermion propagator is not 
in practice subtractively renormalized. Most of 
these studies have assumed an initially massless 
theory and have renormalized at the ultraviolet 
cutoff of the integrations, taking Z\ = Z2 = 1. 
We describe here some results |Q of a study of 
subtractive renormalization in the fermion DSE, 
in quenched strong-coupling QED4. (In the con- 
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text of this study of QED, the term "quenched" 
means that the bare photon propagator is used in 
the fermion self-energy DSE, so that Z3 = 1. Vir- 
tual fermion loops may still be present, however, 
within the vertex corrections.) Results are ob- 
tained for DSE with three different vertices: the 
bare 7 M , the minimal Ball-Chiu vertex form |t]], 
and the Curtis-Pennington vertex (||||lO|,[ll) ■ 

2. DSE and Vertex Ansatze 

The DSE for the renormalized fermion propa- 
gator, in a general covariant gauge, is 

S-\p 2 ) = Z 2 (fi,A)y-m (A)}-iZ 1 (^A)e 2 

x / ^rS(k)r»(k,p)D^( q ); (2) 

here q = k — p is the photon momentum, fi is 
the renormalization point, and A is a regular- 
izing parameter (taken here to be an ultravio- 
let momentum cutoff). We write mo (A) for the 
regularization-parameter dependent bare mass. 
The physical charge is e (as opposed to the bare 
charge eo), and the general form for the photon 
propagator is 

with £ the covariant gauge parameter. Since we 
will work in the quenched approximation and the 
Landau gauge we have e 2 = eg = 47rao and 

D^{ q )^D^{ q )=(-g^ + ^pi^, 

for the photon propagator. 

The requirement of gauge invariance in QED 
leads to the Ward-Takahashi Identities (WTI); 
the WTI for the fermion-photon vertex is 
%THk,p) = S-^k) - S^ip), where q = k-p 
MM. This is a generalization of the original dif- 
ferential Ward identity, which expresses the effect 
of inserting a zero-momentum photon vertex into 
the fermion propagator, dS~ l {p) / dp v — T v {p,p). 
In particular, it guarantees the equality of the 
propagator and vertex renormalization constants, 
Z2 = Z\. The Ward-Takahashi Identity is easily 
shown to be satisfied order-by-order in perturba- 
tion theory and can also be derived nonperturba- 
tively. 



As discussed in this can be thought of as 
just one of a set of six general requirements on 
the vertex: (i) the vertex must satisfy the WTI; 

(ii) it should contain no kinematic singularities; 

(iii) it should transform under charge conjuga- 
tion (C), parity inversion (P), and time rever- 
sal (T) in the same way as the bare vertex, e.g., 
C _1 r M (fc,p)C = — rj(— p, —k) (where the super- 
script T indicates the transpose); (iv) it should 
reduce to the bare vertex in the weak-coupling 
limit; (v) it should ensure multiplicative renor- 
malizability of the DSE in Eq. (|J) ; (vi) the trans- 
verse part of the vertex should be specified to en- 
sure gauge-covariance of the DSE. 

Ball and Chiu have given a description of 
the most general fermion-photon vertex that sat- 
isfies the WTI; it consists of a longitudinally- 
constrained (i.e., "Ball-Chiu") part which is 
a minimal solution of the WTI, and a basis set of 
eight transverse vectors T^(A;,p), which span the 
hyperplane specified by qtl T^{k,p) = 0, q = k—p. 
The minimal longitudinally constrained part of 
the vertex is given by 

T£ c (k,p) = \[A{k 2 ) + A{p 2 )]^ + ^T^x 

|[A(fc 2 ) - A(p 2 )]l±^ [B(k 2 ) - B(p 2 )]} . (3) 

The transverse vectors can be found for example 
in Ref . [[u] . A general vertex is then written as 

8 

T^k,p)=T% c (k,p) + J2nTt{k,p), (4) 
i=i 

where the Ti{k 2 ,p 2 ,q 2 ) are functions which must 
be chosen to give the correct C, P, and T 
invariance properties. Curtis and Pennington 
P|,P|, ^0|JTl]| eliminate four of the transverse vec- 
tors since they are Dirac-even and must gener- 
ate a scalar term. By requiring that the vertex 
r M (fc,p) reduce to the leading log result for k ^s> p 
they are led to eliminate all the transverse basis 
vectors except T£ , with a dynamic coefficient cho- 
sen to make the DSE multiplicatively renormaliz- 
able. This coefficient has the form Te(k 2 ,p 2 , q 2 ) = 
(l/2)L4(fc 2 ) - A{p 2 )]/d{k,p) 1 where d(k,p) is a 
symmetric, singularity-free function of k and p, 
with the limiting behavior lim^-g^ d(k,p) = k 2 . 
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[Here, A(p 2 ) = l/Z(p 2 ) is their l/T{p 2 ).] For 
purely massless QED, they find a suitable form, 
d(k,p) = (k 2 — p 2 ) 2 / (k 2 + p 2 ) . This is generalized 
to the case with a dynamical mass M(p 2 ), to give 
d(k,p) = [{k 2 -p 2 ) 2 + [M 2 {k 2 ) + M 2 (p 2 )] 2 ]/{k 2 + 
p 2 ). 

3. The Subtractive Renormalization 

As discussed in Ref. j| one first determines a 
finite, regularized self-energy, which depends on 
both a regularization parameter and the renor- 
malization point; then one performs a subtraction 
at the renormalization point, in order to define 
the renormalization parameters Z\, Z 2 , Z%. The 
DSE for the renormalized fermion propagator, 

S-\p) = Z a (M,A)[jrf-mo(A)]-E , (M,A;p) 

= i>-m{[h) -E(fi;p) , (5) 

where the (regularized) self-energy is 



A;p) — iZie 2 



d*k 
j2nY 



(6) 



[To avoid confusion we will follow Ref. |5J and in 
this section only we will denote regularized quan- 
tities with a prime and renormalized ones with a 
tilde, e.g. £'(//, A;p) is the regularized self-energy 
depending on both the renormalization point \i 
and regularization parameter A and £(/z;p) is 
the renormalized self-energy.] The self-energies 
are decomposed into Dirac and scalar parts, 
£'(/i,A;p) = Y! d {^A-p 2 )i> + SUM,A;p 2 ) (and 
similarly for the renormalized quantity, E(/i,p)). 
By imposing the renormalization boundary con- 



= i> - m(p), one 



gets 



^, s (M,A;p 2 



dition, S -1 ^) 

the relations 
Y,' d s (n, A; n 2 ) for the self-energy, Z 2 (jU,A) = 
1 + S' d (/i, A; /i 2 ) for the renormalization, and 
m (A) = [m(/i) - S' s (/i, A;/x 2 )] /Z 2 (fi, A) for the 
bare mass. The chiral limit is approached by tak- 
ing m{p) = X' S (^,A;/^ 2 ) = for fixed p.. This 
is the limit in which the renormalized axial vec- 
tor (non-flavor singlet) current is conserved. As 
A — > 00 the numerical results are consistent with 
the bare mass also vanishing as one would expect. 
The mass renormalization constant is given by 



Z m (/j,,A) = TOo(A)/m(/i). The vertex renormal- 
ization, Zi(fj,,A) is identical to Z 2 (fi,A) as long 
as the vertex Ansatz satisfies the Ward Identity; 
this is how it is recovered for multiplication into 
E'(/z,A;p) in Eq. @. Since for the bare vertex 
case there is no way to determine Zi(p, A) inde- 
pendently we will also in this case use Z\ = Z 2 
for the sake of comparison. 

4. Results 

Solutions were obtained for the DSE with the 
Curtis-Pennington and bare vertices, for cou- 
plings ao from 0.1 to 1.75; solutions were also 
obtained for the minimal Ball-Chiu vertex, with 
couplings ao from 0.1 to 0.6 (for larger couplings 
the DSE with this vertex was susceptible to nu- 
merical noise). For the solutions in the subcritical 
range, the renormalization point /x 2 — 100, and 
renormalized masses were either m(/x) = 10 or 30. 
Ultraviolet cutoffs were 1.0 x 10 12 and 1.0 x 10 18 . 
The family of solutions for Landau gauge and 
the Curtis-Pennington vertex with m(fi) = 10 is 
shown in Fig. El Fig. [l] a) shows the finite renor- 
malization Affi); note that for the solutions with 
bare 7", we would have A(p 2 ) = 1 for all values of 
the coupling. Fig. Ej b) shows the masses M(p 2 ). 
Since the equations have no inherent mass-scale, 
the cutoff A, renormalization point /1, m(fi), and 
units of M(p 2 ) or B(p 2 ) all scale multiplicatively, 
and the units are arbitrary. The results for all 
vertices are qualitatively similar. For the Lan- 
dau gauge results we find that as the value of ji 
increases, Z 2 approaches 1, and also that Z 2 is 
almost independent of the cutoff A, at least for 
A»/i. The solutions are extremely stable as the 
cutoff is varied. 

In all cases the mass and finite renormalization 
were stable with respect to very large variations in 
cutoff. The renormalization constants Z 2 (p, A) = 
Zx(fA,A) remain finite and well-behaved with in- 
creasing A in contrast to what happens in pertur- 
bation theory. 

5. Summary and Conclusions 

A full discussion of the Landau gauge results 
can be found in Ref. |l4| . Results have also been 
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p (spacelike) 



Figure 1. Family of DSE solutions with the 
Curtis-Pcnnington vertex and subcritical cou- 
plings a = 0.1 to 0.9 in the Landau gauge. The 
renormalization point is /i 2 = 100, and renormal- 
ized mass is m(/i) = 10 for all cases shown. Ul- 
traviolet cutoffs arc 1 x 10 12 for all cases, (a) the 
finite renormalizations A(p 2 ); (b) the mass func- 
tions M(p 2 ). 



recently obtained in other covariant gauges and 
this work is currently being prepared for publica- 
tion. A careful comparison with lattice formula- 
tions of gauge-fixed QED can now be attempted, 
since we are free to choose the renormalization 
point in this work as the zero-momentum point 
for the n-point Green's functions just as is con- 
ventionally done in lattice gauge theory studies. 
Such comparative studies should yield significant 
benefits. 
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